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Executive Summary 


Multi-body flexible structures, such as those currently under investigation in the spacecraft 
design, are large scale (high-order) dimensional systems. As such, controlling and filtering such 
structures is a computationally complex problem. This issue is particularly important when many 
sensors and actuators are located along the structures that need to be processed in real-time. 

This report summarizes research activity focused on solving the signal processing (that is, 
information processing) issues of multi-body structures. A distributed architecture is developed in 
which single loop processors are employed for local filtering and control. By implementing such a 
philosophy with an embedded controller configuration, a supervising controller ( as simple as a 
workstation or personal computer) may be used to process global data and make global decisions 

as the local devices are processing local information. 

A hardware testbed - a position controller system for a servo motor- is employed to 
illustrate the capabilities of the embedded controller structure. Several filtering and control 
structures which can be modelled as rational functions can be implemented on the system 
developed in this research effort. Thus the results of this study provides a support tool for many 
Control/Structures Interaction (CSI) NASA testbeds such as the Evolutionary model and nine-bay 

truss structure. 
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I. Introduction 


Large space structures are characterized by large order dynamics; typically a subset of these 
flexible modes are selected as significant parameters but even this set can be on the order of 100 
states. Because of the large dimensionality issue, it becomes desirable to process sensor data 
efficiently and effectively in order to send appropriate information to the controller. 

Several researchers have looked at the modelling issues of large space structures [3-7, for 
example]. The classical approach is to identify the significant flexible modes of the system 
experimentally and to assume that little, if any, of the spillover effects will significantly reduce the 
desired closed-loop system response. 

Even with a model reduction of the number of flexible modes, there are usually a large 
number of modes close to the imaging axis that need to be controlled. Compensation requires 
obtaining sensor information, processing the data efficiendy to produce a command signal and 
transmitting the control strategy through the actuator to the system. With the large spatial and 
modal structure of the system to be controlled, many space structures process a large number of 

sensors and actuators. 

Sensing [1-2] and the associated pre-conditioning filter problem [8, for example] are 
important issues that need be addressed efficiently for real-time applications. Typical bandwidth 
requirements for the significant modes of structure are weighed against sensor and actuator 

bandwidth as well as inherent noise problems. 

Towards this end, this research effort has focused on the development of a hardware tool 
with associated software and system architecture that addresses some of the requirements for 
digital signal processing for multi-body flexible structures. The approach is to use a distributed 
architecture whereby fast single loop (sensor-control-actuator) processing can be accomplished 
with a supervisor for global coordination. The overall system structure is now enumerated. 



H. System Architecture 


For a space structure made up of flexible lightweight components, a large amount of data 
need be processed. Classical architectures use the structure as shown in Figure 1. The entire 
structure incl udin g sensors and actuators is handled as one system; all of the sensor data is sent to a 
filter, a state estimator if required and a centralized controller. The controller returns a command 
input to all of the actuators. 

The architecture developed here is based upon a distributed configuration as illustrated in 
Figure 2 and detailed further in Figure 3. Note that in this architecture, the distributed space 
structure is partitioned according to sensor-actuator loops. Each subsystem then requires local 
control around the sensor-actuator pair. There are several ways to parntion a distributed system 
into subsystems. One method spatially breaks up the structure into parts, with interaction 
dynamics between subsystems. Another method partitions the structure according to mode groups 
where sensors, including filters, separate the frequency spectrum of the structure according to 
different bandwidths. 

Given that some partitioning has occurred, the next design step is to develop the local 
controller which is implemented for real-time processing. A central processor (supervisor) is used 
to coordinate global specifications (stability, overall time response characteristics, for example). 

There are several filters that can be employed for pre-conditioning of sensor data. 

Appendix B provides a summary of several designs using the Butterworth, Chebyshev, Elliptic 
(Cauer) and Bessel filters on sensor data from the CSI Evolutionary Model. For typical space 
structure system specifications, the following observation can be made about these filters (for a 
given order): of the structures under study, the Butterworth Filter in the simplest filter to design; 
the Cauer filter gives the best cut-off characteristics but the worst phase distortion (nonlinearity). 
Hence the Butterworth and Cauer filters are recommended for further study as signal processors 
for pre-conditioning. 
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Table 1. 



Filter type 

Passband Characteristics 

Cart-off Space 

Stop band 

1 . 

Butterworth 

Maximally Flat Magnitude 
Moderate Phase Distortion 

Moderate Cut-off 

Moderate 

Attenuation 

2. 

Chebyshev 

Ripples in Passband 
Extreme Phase Distortion 

Sharper Cut-off 

Better Attenuation 
than Butterworth 

3. 

Cauer 

Ripples in passband 
Extreme Phase Distortion 

Sharper Cut-off 
than Chebyshev 

Ripples in 
Stopband 

4. 

Bessel 

Flat Magnitude. Less Phase 
Distortion than Butterworth 

Slow Cut-off 
Specifications 

Less Attenuation 
than Butterworth 


With the filter design selected, an appropriate controller is required for local control. Both 
static dissipative control [12] and virtual passive control [13] have been studied. A discrete version 
of the virtual passive control has also been developed [14]. The controller and filter structure are 
assumed linear in nature and thus can be described by a transfer function representation or ARMA 
(auto-regressive, moving average) architecture. With this in mind, the local filter and controller is 
implemented on an embedded controller, with a higher level computer (minicomputer or 386-based 
computer, for example) as a supervisor (Figure 4). This controller is based upon the Dallas 2250 
chip and associated interface. 

The DS2250 Control system is a microcontroller-based, adaptable controller system that 
can be changed to control a variety of systems. The main idea behind the system is that there is a 
sensor input, a transfer function representing the processor and an actuator output The 
microcontroller can be attached (via a serial cable) to a host computer to allow a variety of transfer 
functions to be tested. 

The basic system contains a DS2250, Intel 8051 -based, microcontroller. This controller 
has 32k of battery backed-up RAM at it’s disposal, as well as four bi-directional, parallel I/O ports 
for use with peripherals. The host computer has a C cross-compiler and a download program that 
works with this micro-controller. 
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Attached to the DS2250 are an RS-232 interface chip, an Analog to Digital converter, and a 
Digital to Analog converter. The RS-232 chip is used for communication with the host computer 
for down- loading and control procedures. The A/D and D/A converters are used to take the signals 
from the sensor, and send the signals to the system under study (Figure 5). 

The software development of this research effort includes several test programs for the 
converters and the serial port, as well as a library of functions that access the converters. Further, 
a front-end program has been developed which allows the user to enter the filter/controller 
characteristics which are then down-loaded to the DS2250 system. Once this is done, control is 
passed on from the supervisor to the local controllers for local signal processing. These functions 
are now discussed. 

HI. The DS2250 Embedded-Controller Chip 

The heart of the distributed processor unit (DPU) is the Dallas DS2250 embedded 
controller (Figure 6). The 40-psi "Soft Micro Stik" is an Intel 8051 -based micro controller which 
has 64K bytes of nonvolatile SRAM for program and data. Using an 8-bit C-MOS baseline, the 
DS2250 allows the user to change program memory at any time. Serial loading uses an on-chip 
I/O part to accept data from a host computer via an RS-232 part. This allows a user-friendly 
program development feature in which the designer can develop and modify code using a personal 
computer and download pertinent information for control to the DS2250. It is this feature that 
makes the DS2250 an attractive device for distributed processing development. 

Other features of the DS2250 include: 

( 1 ) data RAM or data registers retain information in the absence of a voltage reference, 

(2) extensive security features allows the designs to control unauthorized usage of 
program RAM, 

(3) the DS2250 is compatible with several existing software and hardware support for the 
Intel 8051 -embedded controller, 
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(4) the instruction set also includes interrupts, event counters and timers, serial I/O port 
capability for synchronous and asynchronous operations and special function 
registers, 

(5 low power requirements are maintained by power cycling during idle times, 

(6) automatic restart if an errant piece of software was executed; provides the program 
with robustness to power failure or power surges. 

The DS2250 belongs to a family of Dallas embedded controller chips which include the 
DS5000 (a parallel to the DS2250 but with less RAM), one DS5001Fp (an enhanced version of the 
DS5000 but with a 80-pin quad flat package) and the DS2256 (the lowest-power version of the 
DS5001 with 128K bytes of nonvolatile RAM). All of these chips offer flexible read/wnte 
memory, SRAM, which is much more attractive than ROM-based controllers. A detail description 

of this chip now follows. 

The Dallas DS2250 (Figure 7) is a family of 8-bit CMOS microcontrollers [10]. The 
DS2250 is implemented on a 2.65 in by 0.84 in 40-in SIP. The DS2250 is available in various 
memory sizes, up to 64K of CMOS RAM, in two maximum clock speeds, 12MHz or 16MHz, and 
with or without a permanently powered clock/calendar. The RAM memory is nonvolatile, due to 
an onboard Lithium battery that supplies sufficient power to the RAM for memory retention when 
power is removed to the system. The DS2250 uses a 8051 instruction set. The DS2250 requires 
only a crystal and a 5V power supply to become an operational system. 

The DS2250 has numerous features that allow for easy implementation of the system. All 
of the CPU registers are mapped as Special Function Registers (SFR s) and they are identical in 
number and function to those of the 8051 (i.e., the Accumulator, the Stack Pointer, the B Register, 
the Program Status Word and the Data Pointer Register). There are 128 scratchpad registers that 
can be accessed by the CPU. The SFR’s are used to establish several other important features of 
the DS2250. 

There is serial I/O capability in the DS2250. The input and output buffers are accessed as a 

single SFR (SBUF) and control is performed a serial control SFR (SCON). Data that is to be 
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transmitted is loaded into the SBUF, while received data is read from SBUF. Format of the serial 
data, such as the number of data bits, the parity bit and the number of stop bits, is set by the 
control word that was programmed into SCON. The transmit and receive lines for the serial 
interface use the two least significant bits of the 8-bit parallel port, P3. The band rate is 
programmed by using one of the programmable timers. The output of the serial I/O port is at 
CMOS levels (supply voltage is 5V), and thus for compatibility with RS232, these levels need to 
be converted. 

There are two programmable 16-bit timers/counters. Each timer has a SFR for each the 
high and low byte of the 16-bits. These four SFR can be written to set the timer's intervals, or 
read to determine the number of events that occurred related to that particular timer/counter. The 
timers operate by incrementing the timer’s 16-bit register each time a pulse is supplied to the timer, 
providing the timer has been enabled. The pulse can come either from 1/12 of the CPU clock 
oscillator or from off the board, bits 4 and 5 of parallel port P3. Source of the pulse, timer 
enabling and mode of timer operation is set by the timer control SFR (TCON). There are various 
modes in which the timer can operate (i.e., auto-reload, 8 or 16-bit operation) which are set by the 
contents of TCON. TCON also controls the generation of an interrupt by the timer. 

The DS2250 has four 8-bit bidirectional parallel I/O, denoted by PO, PI, P2 and P3. These 
ports also have other associated functions. P3 has been seen to have two lines associated with 
serial I/O and two lines associated with the timers. Two more lines of P3 have additional use for 
external interrupts and the remaining two lines are associated as timing lines for expanded memory. 
If any of these alternate uses for the data lines are not used, then they revert to parallel I/O. Ports 
PO and P2 have the alternate use of forming a 16-bit data bus for expanded memory (off SIP 
memory). If the expanded memory is not implemented PO and P2 remain as parallel I/O ports. All 
of the ports have pull ups on each data line. To output data, the data is written to the desired port. 
While for input, a OxFF is first written to the port and then the port is read. 

There are various configuration parameters that are protected from errant modification. 
Some of these are the memory partition data bits, the security lock bit, the watchdog timer enable 
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and reset bits, the stop mode bit and Power-On-Reset bit. To protect this data, a timed access 
method is used to alter any of these data bits. To modify any of these data bits, a OxAA is written 
to the Timed Access Register followed within two machine cycles of a write of 0x55 to the Timed 
Access Register. The next instruction can alter any of these configuration parameters, within the 
next two machine cycles. If the 0x55 is not written within the required time following the writing 
of OxAA, or the timed access restricted data is not modified within its time limit of the writing of 
0x55, then the modification of the data is not permitted. This method is used to prevent errant code 
from destroying the machine configuration. The logic is that performing the two writes and the 
modification of the configuration in sequence and in the time limit is highly unlikely by errant code 
(i.e., pure chance). Information contain in these configuration parameters can be read without 
going through the timed access routine, but they can not be altered with following the routine 
exactly. 

There is circuitry on the DS2250 to detect a power fail condition. If enabled this will 
generate an interrupt. This will allow for the storing of the CPU state in memory for retrieval 
during a power up. With the power off the system will retain its memory contents because of the 
Lithium battery. This allows for the program to resume operations after the power failure. 

The DS2250 also has a Watchdog Timer. This is used to reset the DS2250 in the event that 
software control is lost The Watchdog Timer is enabled via a Timed Access setting of the 
Watchdog Timer Enable. The timeout period is equal to 122880 machine cycles. The Reset 
Watchdog Timer (RWT) is set high using a Timed Access routine. The RWT must be set high 
before the Watchdog Timer times out, as long as the Watchdog Timer is enabled, otherwise the 
Watchdog Timer generates a reset signal for the DS2250. The Watchdog Timer provides for a 
reset the DS2250 gets into a "Death Grip" with some device, or any situation occurs in which the 
software is unable to reset the timer within the allotted time. 

The DS2250 also has a Security Lock, Encrypting Logic and internal Vector Ram for 
interrupt servicing. The Security Lock provides a lock on the program memory and configuration 
data. The lock is set at the completion of programming, if desired. The program memory can be 
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read only if the Security Lock is reset. Once the Security Lock is set, the program memory can not 
be read and if the Security Lock is then unlocked the contents of the CPU and memory are zeroed. 
The CPU can still operate from the program in memory when the system is locked. At the time of 
programming the DS2250, the program can be encrypted. The encryption routine is based on the 
use of a 40-bit encryption word that is set during programming. The encryption word is used to 
encode both the program code, data and addresses used by the system. The internal Vector RAM 
is used to by the system during encrypted operation to prevent sending of addresses or data to the 
memory chips on the DS2250 SIP during an interrupt service routine or during a power on reset. 

The DS2250 has no ROM for program storage. All program and data memory is in the 
system's RAM, and since the memory is nonvolatile due to the Lithium battery it is not lost even 
when the system’s power is removed. This allows for the changing of the memory contents 
without the need to erase an EPROM or changing to a new DS2250. The memory can be 
programmed either through a serial or parallel interface. During programming, in addition to the 
loading of the program code, the system's memory partition is set, the system can be locked and 
the encryption key word can be set. Parallel programming of the DS2250 has not been perused, 
due to the simplicity in hardware and software available for serial programming. 

IV. The Software Package 

There is a programming package available from Dallas Semiconductor, called K3T5K [10]. 

KIT5K runs on an IBM PC/AT compatible machine operating under DOS2.0 or higher. This 

software communicates via either of the COM ports on the machine to the DS2250. Since the 

COM port is RS232 and the DS2250 requires CMOS logic levels, a RS232 to CMOS convener is 

required. The KIT5K software toggles the DTR line of the COM port which can be connected to 

the Reset and PSEN* line of the DS2250 to put the system in the serial programming mode, or this 

can be performed manually (if desired). The DTR must also be converted from RS232 to CMOS 

levels. Once the Reset and PSEN* lines are activated to assert serial programming mode, the 

DS2250 senses the serial port receiver line to detect the <CR> transmitted by the programming 
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computer. From this serially transmitted character, the DS2250 is able to determine the baud rate at 
which it is to be programmed and the serial programming mode is established. Not all baud rates 
are available. The usable standard baud rates for various crystal frequencies are published in the 
User's Manual DS5000. 

There are numerous commands for the KIT5K software, and all of them will not be 
examined here. Most of the commands can be placed into three categories: PC Control, 
Programming Control and Memory Editing. Some of the commands from each category will be 

examined. 

The PC Control category contains commands primarily used to control the PC that KIT5K 
is operating on. The COM command directs through which port the KTT5K software will 
communicate with the DS2250. The SPEED command sets the baud rate and DTR toggles the 
DTR line, of the selected com port Since on each entry into KTT5K requires that the com port and 
the baud rate to be set, to avoid exiting KTT5K there is a DOS command that allows DOS 
command execution without leaving KTT5K. In addition, there are some DOS commands built 
into the KTT5K to allow changing the current directory of listing its contents. 

The Programming Control commands are to facilitate programming of the DS2250. 
Obviously, the dominant command is PROGRAM which performs all the steps required to 
program the DS2250 with a specified Intel hex format file containing the program code. There are 
individual commands to perform each item required to program the DS2250 from the PC. The 
PROGRAM command automatically performs the required items and queries the user about 
programming options. Once KIT5K has been initialized, PROGRAM along with the specified hex 
file is used as a command. KIT5K then executes a DTR command to initialize serial programming 
and then a <CR> is transmitted to establish serial programming. An UNLOCK command is 
executed to unlock the DS2250 Security Lock. The software then queries the user for 
configuration data, namely memory size, partition address, encryption and key if needed, 
beginning address, ending address and if the programming is to be verified. After these values 
have been entered, the PROGRAM command then executes the various commands needed to 
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program the DS2250 as specified. The RANGE command is used to set the memory size of the 
DS2250. The PARTITION command is executed to set the specified address for the program/data 
space partition. If encryption is being used the KEY command is executed to load the DS2250 
with the encryption key. The hex format file is then downloaded to the DS2250 via the LOAD 
command. If verification of the downloaded program is required the VERIFY command is 
executed to compare the hex file with the memory contents. If the system is to be locked, a LOCK 
command is executed to assert the DS2250's Security Lock. And finally, the RUN command is 
executed which places the DS2250 in the execution mode. In order to simplify the process a batch 
and configuration file can be used to perform these tasks. The batch file delivers the commands to 
be executed, while the configuration file has the system configuration data in it. 

The Memory Editing commands allow for the altering and examination of the DS2250 s 
memory. These commands are similar to the commands of DDT or DEBUG. The FILL command 
fills the specified memory locations with the specified constant The EDIT command permits the 
user to sequentially view and alter memory starting at the specified address. DUMP stores the 
contents of the DS2250 memory between the two specified address in a specified file. The 
UNASM command unassembles the contents of the DS2250's memory. 

There are numerous ways of generating an Intel hex format file containing the program 
code. There are many assemblers and cross compilers that run under DOS that produce these files. 
The Franklin 'C' Compiler [11] does meet this requirements and operates under DOS3.0 or higher. 
This 'C* compiler and its associated programs generate Intel hex format files for the 8051 family of 
microcontrollers and related microcontrollers. The compiler package is composed of the C 
compiler, an assembler, a linker, a library manager, and an object to hex converter. 

The 'C compiler implements an ANSI standard 'C'. There are extensions to this package 

to allow for easy access to the microcontroller's hardware. This is accomplished to two ways. 

First, the specific Special Function Registers are defined in a CPU specific header file that contains 

the names of the SFR/s mnemonic as a variable. Thus, an output of 0xF8 by parallel port 2 is 

achieved by the statement, P2=OxF8;. One of the extensions of the compiler allows for bit 
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communications with the SFR's mnemonics followed by a " i, with i being the desired bit number. 
To reset the third bit of port 2, use, P2 A 3=0;. Operations with the SFR’s or particular bits of the 
SFR's is easily performed in the 'C program. In addition particular bits can be given unique 
names by the define statement. The 'C program can access all of the DS2250, without the user 
reverting to assembler code. It should be noted that a timed access write and interrupt servicing 
can be performed in 'C' code. 

The Franklin 'C compiler, while it can easily communicate with the DS2250's hardware, 
does not permit for the inclusion of assembler code inside the 'C code. Items that need to be 
executed in assembler must be written as an assembler function, assembled and linked to the object 
code generated by the compiler. In addition, the compiler can be linked to functions generated m 
Intel PL/M 51. The output of the tinker has a file extension of abs. This absolute object file may 
contain symbol information for symbolic debugging. The object-hex-symbol converter generates 
an Intel hex format file from the output of the tinker. This hex file is then programmed into the 
DS2250 by the KIT5K software. 

Besides these system files, a batch file has been developed as a user-friendly 
communication between the designer and the DS2250 system. This file, tf.bat, is resident on the 
supervisor system (IBM system in this study) and allows the user to create a continuous-time 
transfer function or a discrete-time (digital) transfer function. In the continuous-time case, a 
discrete-time equivalent transfer function is created using the bilinear z-transformation. Currently, 
the system allows for a fifth-order continuous-time system or a fifteen-order discrete-time system. 
These parameters can easily be changed by an appropriate modification in the driver software. 

Once the local processor coefficients are entered through the supervisor system, it is down- 
loaded to the DS2250 embedded controller. From this point on, the local controller takes over. 
Through the use of data acquisition chips, the local processor obtains information from the sensor 

and sends commanded inputs to the actuator. 

Note that the transfer function form has applicability to many processing systems. The 

filters designed in Appendix B, for example, can be readily tested using this local processor. 
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Further, classical PID (proportional, integral, derivative) control, phase compensation and virtual 

passive control can be implemented on this system. 

Issues such as bandwidth limitations, system order, voltage level constraints and software 

requirements are presented in Appendix C. 

V. Data Acquisition Hardware 

The analog-to-digital converter is the Analog Devices ADCO804C chip. This system has 
8-bit resolution with differential inputs and a sampling rate of up to lOKHz. The converter has a 
reference voltage of 5v and operates with a 256-resistor DAC network. 

The digital-to-analog converter is the Analog Devices AD7530/533. This system is a 10-bit 
multiplying DAC using CMOS technology and uses a 5v reference (it can operate with voltage 
references up to 15 volts). An inverted R-2R ladder is used for the binary weighed operations; a 
sampling rate is typically quoted at 160ns (6.7MHz). The DAC used in this study is set for bipolar 
operation (4-quadrant multiplication). 

Finally, an RS-232 interface circuit (the ICL232CPE) is used to interface the Dallas 
DS2250 controller to an IBM computer. Pin connections are shown in Figure 8. All of the data 
acquisition components are fully functional and are transparent to the user. Further software to 
analyze the data acquisition units and to transfer data and command signals back and forth between 
the embedded controller and the data acquisition system have been fully tested; upon start up, these 
modules are activated and the user need not access any additional files. 

A discussion of the start up procedure and hardware issues is provided in Appendix C. 

VI. Example-Using the DPU for Motor Control 

A two-beam flexible testbed has been constructed for collocated and non-collocated control 
studies. A separate report [15] summarizes some of the initial studies on non-collocated control 
using zero placement and phase compensation. The testbed has a stepper motor produced by 
Servo Systems which rotates along the z-axis (vertical). Two flexible steel alloy beams are 
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attached to a mounting unit These beams have been selected so that all movement is in the plane 
parallel to the ground; little torsional motion is generated, as desired. 

Two servo motors (produced by the Maxon Company) provide actuation for these beams. 
These motors were selected because of the speed requirements to suppress vibration. Two Copley 
control servo amplifiers provide the necessary voltage from the computer to the motor. Each motor 
has an incremental optical encoder (produced by Hewlett-Packard) with 1024 cycles per resolution. 

Each beam has a 350 ohm strain gauge mounted along a side which sends non-collocated 
position information through an amplifier circuit to the data acquisition system. Figure 9 illustrates 
the two-beam flexible testbed. 

For the purpose of demonstrating the DPU, the two-beam unit was replaced by a single 
beam with a different mounting unit. The objective of the test was to select controller gains to 
suppress beam vibration. Proportional-derivative control was selected for illustration. Files 
download.bat 9.1, 9.2, 9.3 contain the data for these trials. File 9.1 contains the proportional 
gain, 9.2 contains the derivative gain and 9.3 contains the desired terminal angle in radians. 

First, the user loads the start-up file tf.bat. Then the software prompts the user for 
controller/filter parameters and upon completion of the data entry, downloads the data and passes 
control to the DPU. 

Figure 10 illustrates the results for four different runs. For this single study, one observes 
that a proportional gain of -10 and a derivative gain of +8 provides a desirable result for vibration 
suppression. The plots are in seconds and show the output of the strain-gage, i.e., the position of 
the beam during the movement from 0° to 45*. 

Other tests have been performed on the DPU and it is noted that this embedded controller 
system is a viable architecture for distributed signal processing application for flexible structures 
such as those currently being considered in spacecraft designs. 

Bandwidth requirements for most applications under study are satisfied with the lOKHz 
data acquisition unit (a higher bandwidth can be easily added on by changing the DAC and 
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adjusting the clock for the A/D); further adequate accuracy is achieved into the 8- 10 bit data 
conversion and 8-bit microcontroller. 

VII. Concluding Remarks and Areas of Future Work 

This research effort has resulted in an embedded controller system which can be used as 
part of a distributed signal processing architecture for flexible multi-body structures currently under 
investigation at NASA. The distributed processing unit (DPU) is a self-contained processor which 
can be attached to a supervisory computer through an RS-232 port and to a sensor-actuator pah- 
through a data acquisition system. Accuracy and sampling rates are satisfied for many multi-body 
structures currently being used at the Spacecraft Dynamics Branch. 

Areas of future activities are many. First of all, microcontroller chips are constantly being 
upgraded and many of the new designs are smaller and require less power. For example, the 
Dallas DS2255 which is currently under design/fabrication is a single 68-Pin SIP stik data 
acquisition system (A/D and D/A). This chip has 12-bit A/D accuracy, 8-channel input, has a 
digitally controlled channel selector and a 50KHz maximum throughput rate. There are 8 analog 
outputs and an 8-bit output resolution. Further this chip has an addressable serial interface, 
requires just a +5 reference voltage and an external A/D clock. This chip would reduce the size of 

the DPU. 

The DS2255 instrument stik could be interfaced with the DS2256 power miser micro stik - 
an embedded controller chip also currently under design and fabrication. This 68-Pin SIP stick is 
an extremely low power chip (+3v to +5v reference) which includes 128K bytes of nonvolatile 
SCRAM, a watchdog timer, a permanently powered timekeeper, has 29 bits of user-defined I/O 
ports and a kickstarter if an interrupt or idle status was previously enabled. 

Once these chips are available, the existing DPU developed in this study can be directly 
transferred to upgrade the hardware technology. Speed and com compactness would thus be 
upgraded. 
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Another area of future work is in the area of filter designs. Classical algorithms have been 
employed for the pre-conditioning processors. Thus a transfer function model (or ARMA 
structure) has been implemented in the user-friendly driver package. Other filter designs such as 
bandpass, notch and all-pass structures as well as adaptive, self-tunning systems may be 
developed. This would enhance the capabilities of the DPU when the model of the structure under 
study is unknown or certain parameters are uncertain. Further the filter design can be enhanced to 
work with the controller under study. To do this, many controller structures allow some flexibility 
in zero-placement or eigenvector selection. The filter component not only has the requirement to 
suppress certain frequency components but also has the disadvantage of adding plane distortion or 
time delays. The controller has to compensate for some of these problems. 

Finally, distributed processing systems requires some protocols between the supervisor 
and local controllers. Approaches from computer sciences (queuing, token nng structures) and 
optimal control theory (min-max performance measures, game theory) may be employed to handle 
priorities of data sharing between local controllers when necessary. 

These areas of future research would benefit the current objectives of the control/structure 
interaction program by providing an efficient control architecture for multi-body flexible structures. 
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Figure 1 : Overall system with estimator and control blocks 
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Figure 2: Distributed processor architecture 
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Figure 3: A large-scale distributed system 
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Figure 5: The hardware architecture 
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Figure 6: The Dallas 2250 hardware design 
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Figure 8: Bottom view of RS-232 interface components 
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Appendix A 

Driver Program for Processor Coefficients 
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a 


A: \ > 


10 
.20 
30 
4 0 
50 
50 
70 
R0 
^0 
100 
l LO 
120 
1 30 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230 
240 
250 

255 

256 
260 
270 
280 
290 
300 
310 
320 
330 
340 
350 
360 
370 
380 
390 
400 
410 
420 
430 
440 
450 
4 60 
470 
480 
490 
800 
510 
520 
530 
540 
550 


REM 

REM Andrew Meyers 

REM 
REM 
REM 
REM 
REM 
REM 
REM 
REM — 

CLS 


Fall 1991 

This program gets the coefficients either of a cont on 
time transfer function or a d iscr ete-time transfer! Fun- ction 
from a user. If a continuous-time transfer function is entereJ 
it- is converted to discrete-time. Whem ail of lihis is comp 
the function is downloaded to the micr ocontr oler and the mot. 

sequence is begun . 


PRINT 
PRINT 
INPUT 
IF CD 


A 
2 1 


A 2 ONLY" 


key OFF h 

LOCATE 8 , 20 : PR INT " MOTOR CONTROL TRANSFER FUNCTION- 

PRINT : PRINT : PR INT : PRINT : PRINT : PRINT 

PRINT: PRINT 

FOR CONTINUOUS- TIME INPUT 
FOR DISCRETE-TIME INPUT A 
CD 

-1 OR CD =2 THEN GOTO 220 
PR IN T : PRINT "PLEASE INPUT A 1 OR 
GOTO 150 

DIM A( 10 1 ,B( 10 ) , AZ( 30 ) ,8Z( 30 ) 

FOR 1=0 TO 10 
A( I )=0 
B( I 1 = 0 
NEXT I 

FOR 1=0 TO 30 
AZ( I ) = 0 
BZ( I )=0 
NEXT I 

IF CD=2 THEN 1410 

REM 

CLS 

PRINT "CONTINUOUS-TIME 
PRINT : PRINT 
INPUT "Input degree of 
IF ND<=5 THEN 380 
PRINT : PRINT " The order 
GOTO 330 

IF ABS( INT( ND ) ) < >ND THEN 330 
PRINT : PRINT 

PRINT "Numerator coefficients 
PRINT 

PRINT "A(n) + s~n + A(n-l)*s ( n-1 ) + 

PRINT : PRINT 
FOR I=ND TO 0 STEP -1 
PR INT " I nput coefficient A( " 5 I 5 
INPUT A( I ) 

NEXT I 
PRINT : PRINT 
INPUT "Input 


TRANSFER FUNCTION" 
Numerator Poly nomia 1 " ; ND 
of the polynomial must be 


<= 5 " 


are assumed to be of the form: 
. . +A( 1 )*s + A( Oi- 
l- 


Polynomial " : DD 


the degree of Denominator 

IF DD < =5 THEN 530 . . , . ... 

PRINT : PRINT " The order of the polynomial must L. 

GOTO 480 

TF ABS( INT( DD ) )< >DD THEN 480 

PR TNT : PR INT , f form: 

PRINT "Denominator coefficients are assumed to be of the 


160 

-:70 

180 

S90 

5 00 

ilO 

-j20 

d30 

-.40 

:>50 

d60 

= 70 

i80 

o90 

700 

710 

720 

730 

740 

750 

760 

761 
770 
780 
790 
800 
810 
820 
830 
840 
850 
860 
870 
8 80 

890 

900 


PRINT 

PRINT ”B( n)*s'"n + B( n-1 )*s 
PR INT -.PRINT 


( n-1 )+ . . .+B( 1 )*s + 6( 0)" 


O STEP -1 
coef f icient 


B( 


)" 


FOR I -DD TO 
PRINT "Input 
INPUT 8( I ) 

NEXT I 

REN 

CL S 

PRINT :PRINT "DIGITAL CONVERSION WILL NOW TAKE 
PRINT "Input the sampling period (sec) ot the 

INPUT T 
C-2/T 
IF DD=0 
DD=1 
DD=2 
DD = 3 
DD=4 
DO = 5 


PLACE" 
system “ 


IF 

IF 

IF 

IF 

IF 


760 

790 

870 

970 

1090 

1230 


910 
920 
9 30 
940 
950 
960 
9 70 
980 
990 
1000 
1010 
1020 
1030 
1040 
1050 
1060 
1070 
1080 
1090 
1100 
1110 
1 120 
1130 
1140 


THEN 
THEN 
THEN 
THEN 
THEN 
THEN 

REM 

TF-A( O )/B( O ) 

IF DD=0 THEN AZ( 0 )=TF 
GOTO 1380 

REM 7 - 7 

REM First Order bilinear 

REM 

AZZ=B( O ) + B( 1 )*C 
AZ( 0 M A( 0 ) + A( 1 )*C )/ AZZ 
AZ( 1 )=( A( 0 )-A( 1 )*C )/AZZ 
BZ( 1 )-( B( 0 )-B( 1 )*C )/AZZ 
GOTO 1380 

REM 777 

REM Second Order bilinear 
REM 

AZZ=B( O ) + B( 1 )*C+B( 2 )*C 2 

AZ( 0 X A( O )+A( 1 )*C + A( 2 X~2 )/AZZ 

AZ( 1 )=( 2*A( 0 )-2*A( 2 )*C~2 )/AZZ 

AZ ( 2 H A( 


transformation 


transformation 


BZ( 1 )=( 2*bCo )-2*B( 2 )*C 2 )/AZZ 
BZ( 2 )-( B( 0 )-B( 1 ) + C + B( 2 )*C 2 )/AZZ 

GOTO 1380 

REM 

REM Third order bilinear trasformation 

AZZ-B( O ) + B( 1 )*C + B( 2 )*C 2 + B( 3 )*C 3 ^ 

' Z( 0 )-( A( 0 ) + A( 1 )*C+A( 2 )*C~2 + A( 3 )*C 3 )/AZZ 

AZ 1 W 3*A( 0 )+A( 1 )*C-A( 2 )*C~2~3*A( 3 )*C~3 /AZZ 

AZ( 2 ) = ( 3*A( 0 )-A( 1 )*C-A( 2 )*C~2 r3*A( 3 )*C 3 )/AZZ 

AZ( 3 W A( 0 )-A( 1 )*C + A( 2 )*C~2-A( 3 )*C~3 )/AZZ 

BZ i )=( 3*8( 0 )+B( 1 )*C“B( 2 )*C~2-3*B( 3 )*CT3 /AZZ 

BZ( 2 )=( 3*B( 0 )-B( 1 )*C-B( 2 )+C~2+3*B( 3 )*C 3 )/AZZ 

BZ( 3 ) = ( B( 0 )-B( 1 )*C + B( 2 )*C"2-B( 3 )*C 3 )/AZZ 

GOTO 1380 

REM Fourth Order bilinear transformation 

AZZ-B( 0 )+B( 1 )*C + B( 2 )*C~2 + B( 3 )*C 3 + B( 4 )*C >. 4 — . VAZZ 
AZl O ) = ( A( 0 )+A( 1 )*C + A( 2 )+C~2 HA( 3 )*C 3+A( 4 )+C 4 )/AZZ 

AZ( 1 X 4*A( 0 )+2*A( 1 )*C-2*A( 3 4 >* C 4 V 
AZ( 2 )-( 6*A( 0 )-2*A( 2 )*0 2+6*A( 4 )*C 4 )/AZZ 


1210 

1220 

1230 

1240 


' ./• a )/AZZ 

?70 BZ l =(“Wo)t2«B(l)*C-2«B(3)rt'3-4»B(4K 4 )/AZZ 

1 US Sz( 5 8 US S ? & ^34^4 )*C-4 VAZZ 

1200 BZ(4)=(B(0)-B(1)*C + B(2)-KC^2~B(3)*C 3 + B( 4 )*C 4 )/AZZ 

GOTO 1380 

REfl Fifth Order bilinear transfer ma nt ion 

Izso AZZ=B( 0 ) + B( 1 )*C+B( 2 ).C* ‘2 + « 3 )»C :*3 S^C'S )/AZZ 

1260 AZ( 0 )=( A( 0 )tA( 1 )*C+A( 2 )*C *: ' , 3 w c - 3 _ 3 »a( 4 )*C~4-5*A( 5 )*C~5 )/AZZ 

1270 AZ( 1 )-( 5 *A( 0 ) + 3*A( 1 )*CFA (2 )»C 2-A 3 )*C 3 3»A . * > 5 )*c~ 

1280 AZ< 2 )»( 10 *A( O )*2*A 1 *C-2«A 2 W ? »C'3 + 2*A 4 *C"4- 10*A( 5 )»C' 

1 290 AZ( 3 )=( 10*A( 0 )-2»A 1 >*<--2*A ( 2 )*C .t2»A(_3 )*C i > : 2 1 > ( 5 )»C'5 )/AZZ 

1300 AZ( 4 )-< 5 «A( 0 )-3*A( 1 )*C+ A( 2 )*C 2t A( 3 )*C 3 3 « ^ £ j, c - s )/AZZ 
1310 AZ( 5 ) = ( A( 0 )-A( 1 )*C + A( 2 )*C 2"A(3 ) C ( * 4 ^^g*^ 5 )* C ~5 )/AZZ 

1320 BZ( 1 M 5 *B( O ) + 3*B( 1 )*C + B( 2 )*0 “ ( 3 )* c ~ 3 +2 *B( 4 )*C~4 + 10*B( 5 )*C~ 

iSS ss -3 a ius 0 r-ss : g; < 

REM 

REM Transformation complete 

GOTO 1680 

REM Enter the digital transfer function 

PRINT "DISCRETE-TIME TRANSFER FUNCTION" 

440 PRINT-PRINT a 1 11 *Mn 

450 INPUT "Input degree of Numerator polynomial , ND 

IF ABS( INT( ND ) )< >ND THEN 1440 

PRINT : PRINT • • 1 re acenmAfl to be in the form: 

PRINT "Numerator coefficients are 

PRINT , . ^ . 

PRINT " a( O )+a( 1 )*z~-l + a( 2 )*z -+ • 

PRINT : PRINT 

FOR 1=0 TO ND „ . 

PRINT" Input coefficient a( . I . ) 

INPUT AZ( I ) 

NEXT I 
PRINT : PRINT 
INPUT " I nput 
TF ABS( INT( DD ) )<>DD 
PRINT :PRINT . 

PRINT "Denominator coefficients are 

pmr'nM 1 )*z--itb ( 2 )«--2t . . .tW K-i »z--( K-! )tb( k)«‘-k- 

PRINT : PRINT 
FOR 1=1 TO DD 

PRINT "I nput coefficient b( " ; I ; ) 

INPUT BZ( I ) 

NEX T I 

micro- 

REM function. We must now send these coefficient- 

REM controler inorder to dictate the respose ot the motor. 

CLS 

PRINT “The^digital transfer function coefficients: 


5 )/AZZ 
5 )/AZ Z 


5 )/AZZ 
5 )/AZZ 


L 370 
L 380 
L 390 
L 400 
L 4 1 0 
L 420 
L 430 


.460 
L 470 
. 480 
L 490 
L 500 
L 5 10 
L 520 
l 530 
1540 
1550 
1500 
1570 
1580 
1590 
1 600 
L 6 10 
1620 
1630 
1640 
1650 
1660 
1670 
1680 
1690 
1700 
1710 
1720 
1730 
1740 


assumed to be it 
. . +a( k-1 )*z 


( k-1 ) + a( k )*z~- k 1 


t he 


degree of t he 
THEN 1560 


Denominator Polynomial DD 


assumed to be in the form: 


; az( i ) 


1750 

1760 

1770 

1780 

1790 

1800 

1810 

1820 

1821 

1830 

1850 

1860 

1870 

1880 

1890 

1900 

1910 

1920 

1930 

1940 

1950 

1951 

1952 

1953 
1960 
Ok 


print 

FOR 1^ 

PRINT “ 

NEXT I 

PRINT 

FOR 1 = 

PRINT “ 

NEXT I 

INPUT " 

OPEN " 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

REM 

REM 


0 TO 5 

<a( " ; i ; " ) = " 


1 TO 5 
b( " ; i ; ' 


) = " ; BZ( I ) 


enter the value of theta zero 
coeff.h" FOR OUTPUT AS HI 
Hi 


; THZ 


t n 
hi 
m 

«i 

HI 

HI 

HI 

HI 

HI 

HI 

HI 

add 

o f 


AZ( 0 ) 

AZ( 1 ) 

AZ( 2 ) 

AZ( 3 ) 

AZ( 4 ) 

AZ( 5 ) 

BZ( 1 ) 

BZ( 2 ) 

BZ( 3 ) 

BZ( 4 ) 

BZ( 5 ) 

" ; THZ 

more define statements here inorder 
the digital transf er function . 


, " Hdef ine 

aO ” ; 

, "Hdef ine 

al " ; 

, " Hdef i ne 

a2 M ; 

, "Hdef ine 

a3 " ; 

, “ Hdef i ne 

a4 “ ; 

, “ Hdef ine 

a5 " ; 

, " 4 def ine 

bi “ ; 

, " Hdef ine 

b2 “ ; 

, "Hdef ine 

b3 " ; 

, "Hdef ine 

b4 “ ; 

, "Hdef ine 

b5 " ; 

, " Hdef i ne 

thzer o 


SYSTEM 


to raise the order 


type da ta i o . h 

« dataio.h: This include file is designed to make it easie 

* to read and write analog data 


init serial() /* This initializes the proper regist 
/* for using the builtin serial port. 
SCON = 0x52; /* Serial control register */ 

TMOO = 0x20; /* Timer mode type registei 

TCON = 0x68; /* Timer control register */ 

TH1 - 0xf3; /* Timer high byte register */ 


er s */ 
*/ 


. ni t_ad( ) 


unsigned char _por 
PI ]= OxOE ; 

_por tl = PI & Oxfl; 
Pi - _portl ! 0x04; 
PI = _portl ! 0x0c; 


/* This starts up the AD converter by 
/* pulling /wr low for a time, and by 
srtl; /* grabbing and holding /cs low * 


* / 
*/ 

/ 


/% /cs /wr /rd are high */ 
/* pull /cs and /wr low */ 
/* let go of /wr %/ 


isigned char read_ad( ) 


unsigned char _portO, _portl 
_por tl = PI & Oxfl; 

PO = Oxff; 

Pi = _portl ! 0x08; 

_por to = PO ; 

PI = _por tl [ 0x0c ; 
r etur n( _portO ) ; 


/* This function reads a byte */ 
/* from the internal latches of */ 
/* the AD converter . * / 


/% Pull /cs and /rd low */ 
/* Read port 0 
/* let go of /rd */ 


ite_da( _out ) 
signed int _out ; 

unsigned char _portl , 
_put &= 0x03 FF; 

_portl = ( ( _out 
_por t 2 - _out & 

PI = _por t 1 ; 

P2 = _por t2 ; 


/* This function 
/* integer to the 


write a 10-bit * / 
DA converter */ 


_por t 2 ; 


& 0x0300) >> 4) 

OxFF; /* We set 
/* bits 4&5 
/* latch in 


i ( PI & OxCF ) ; 

up the MS bits of the integer 
of port 1. We really need a 
between the ports and the DA 


on 


*/ 


10-bit */ 
co nv . * / 


A : \> 


c 51 COMPILER V3.07, 
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09/11/91 00:00:20 PAGE 1 


DOS C51 COMPILER V3.07 
OBJECT MODULE PLACED I 
COMPILER INVOKED BY: C 


COMPILATION OF MODULE TFCNTRL 
N TFCNTRL .OBJ 

: \C51\BIN\C51 .EXE TFCNTRL .C 


stmt 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 


1 eve 1 


1 

1 

0 

2 

1 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

2 

0 

t— 

y 

1 
1 


source 


^pragma SMALL 
^pragma PRINT 
^pragma CODE 

# include <reg5000 .h> 
# include <stdlib.h> 
^include <stdio.h> 
{♦include <dataio.h> 
{♦include <coeff .h> 


/* register definitions for 8051 CPU 
/* Has atof , and other functions 
/* stardard i/o definitions 
/* a/d and d/a setup 

/* control coefficients from trans4 


*/ 

*/ 

*/ 

*/ 

*/ 


{♦define FALSE 
{♦define TRUE 
{♦define DELAY 


0 

1 


60 


delay! ms ) /» this routine is used to set up delay loops »/ 

unsigned int ms; 

{ .... 

unsigned int i , J > 

for( i=0; i <ms ; i++ ) ( 

f or( j=o ; j < delay ; j++ ) 5 
) 

} 

/ ,„«„„,««*«»»«.*****«***"«**""«*“ n **”**” ,I ‘““*"/ 

/ | lepr Tnouted Transfer-Function Control 

£***************************^ 

/# The controler is placed inside the while loop 

z The a/d reads from 0 to 255 decimal 


*/ 

*/ 

*/ 


ma i n( ) ( 


unsigned char portO; 
unsigned int out , i ; 

float th ,v ,vo ,aa ,bb ,y [10] ,u [10] 

PI = 0x6e; /* PI 

i ni t_ser ia 1( ) ; 
init_ad( ) ; 


putchar( ’O’); 


for ( i = l ; i < =6 ; i t+ ){ 
u [i] =0 . ; 

y [i] =° - > 

> 


/* initial/ set previous values to zero */ 


while ( 1 ) { 



0) ro 


2 

54 2 

55 2 


delay( 67 ) ; 

por tO=read_ad( ); /* read a value from the a/d */ 
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( ( float )por to* .0246 )-( float )thzero ) , 

for ( i =2 ; i <=8; i++ ) 

u[i]=u[i”l] ; 

u [l] =th ; 

/* IF HIGHER ORDER DIGITAL TRANSFER FUNCTIONS RRE REQUIRED ./ 
/* SIMPLY ADO ON TO THE END OF aa ANU 

aa-aO*u [1] + al*u [2] + a2*u [3] + a3*u [4] +a4*u [5] + a5*u [6] ; 
bb=bl*y [l] +b2*y [2] +b3*y [3] +b4*y [4] +b5*y [5] , 
v=aa-bb » 

if ( v < -10. ) 

vo-*" 10 . > 

else if ( v > 10.) 

vo=10 . ; 

else 

vo=v ; 

out = ( int )(46.05*( 10.-vo)); 
w rite_da( out ); /» urite an integer to the d/a «/ 

for ( i«2; i <=S; i++ ) '* update previous values */ 

y[i]=y[i-l] ! 

y [l] =vo ; 
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ABSTRACT 

Multi-body flexible structures such as those currendy under investigation in spacecraft 
designs, are large scale dimensional systems. As such, controlling and filtering issues for such 
structures are computationally complex in nature. One important issue of concern is the effects of 
filtering or pre-conditioning of the sensor data prior to application of the controller strategy. Phase 
distortion, resulting in nonlinear time delays, and spillover effects, due to finite order model and 
control approximations, may lead to undesirable results. This paper investigates several classical 
low pass filter designs under a new light, that is, with its effects on space structure control. Two 
control strategies - static dissipation control and virtual passive control - are investigated and the 
filter/controller configuration is applied to a numerical model of the control/structures Interaction 
Test Article - The Evolutionary Model - which is an experimental hardware currently being studied 
at NASA Langley Research Center. 


I. Introduction 

The next generation of spacecraft will be larger than current spacecraft and will require 
more sophisticated design, construction and operating methods. A sketch of a possible 21st 
century spacecraft is shown in Figure 1. Characteristics common to this class of spacecraft are 
long trass-type bodies with several appendages such as large solar arrays, radiators and circular 
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Figure 1 : Mission to planet Earth geostationary platform 
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antennas. This type of spacecraft will probably contain a large number of lightly damped and 
closely spaced flexible modes at low frequencies. Traditional rigid-body control of spacecraft will 
not be adequate for these spacecraft. Rigid-body maneuvers such as attitude adjustment and 
slewing of appendages may excite flexible modes causing degraded performance and even 
instability or damage to the structure. This has motivated the study of active control methods for 
vibration suppression of large flexible space structures (LFSS). 

The National Aeronautics and Space Administration (NASA) has developed the Controls 
Structure Interaction (CSI) program to address the problems associated with the design and control 
of LFSS. At the Langley Research Center's Spacecraft Dynamics Research Lab (SDRL), a testbed 
has been developed for system identification and closed-loop control experiments of CSI models. 

One challenge regarding the control of LFSS is that accurate numerical models of large 
complicated structures are difficult to obtain. The controls designer often develops control laws 
which incorporates techniques for dealing with errors in the design model, and noise and 
disturbances in the input signal. Most of these techniques require some knowledge of the 
ma gn itude and bandwidth of these uncertainties. However, the complex, multibodied structures 
planned for future missions will make it increasingly difficult to predict these uncertainties. Hence 
it would be ideal to use filters to condition the sensor output so that control design would be less 
demanding. This is the motivation for a study of signal processing techniques, more specifically, 

low-pass sensor filters, for closed-loop control of LFSS. 

The technology for conventional low-pass filters is well developed; they are well 
characterized, and easy to implement With system identification and closed-loop control 
experiments, a first-order sensor filter is often used. In initial attempts to design a controller for 
the CSI Evolutionary model (CEM), a first-order low-pass Bessel filter was used to attenuate the 
response of a 34 Hz mode which was going unstable. The use of this filter allowed for a 40 
percent increase in feedback gains before causing instability. Digital filtering was used in system 
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identification of another CSI test hardware, the Mini-Mast [1]. In both of these cases, the 
frequency above which attenuation is desired, or cutoff frequency, is set high so that there is little 
magnitude or phase distortion in the frequency range of interest 

Generally, higher-order filters are not used in closed-loop applications because the phase 
distortion negatively affects system stability. One example is an attempt to use multiple bandpass 
filters to isolate the modes to be controlled [2]. The frequencies of the modes to be controlled were 
close enough so that the rolloff from one bandpass filter affected the neighboring modes. This 
coupling between the filters and the modes produced system instability. 

With the use of digital filtering, modifications to conventional filters can easily be made, 
and the use of higher-order filters may be possible. In research conducted on sensitivity of the 
Space Station Freedom’s Alpha Joint controller to variations in the structural modal parameters, a 
modified fourth order Butterworth filter is used to attenuate the high frequency structural modes 
[3]. The modification consists of the addition of two zeros in the numerator of the filter transfer 
function to offset the phase distortion caused by the conventional filter. The filter dynamics is 
included within the controller optimization and the filter cutoff frequency is a design variable. 

n. Controls Structures Interaction (CSI) Program Test Article 

The test article selected as the reference problem for this paper is the CSI Evolutionary 
Model (CEM). The CEM testbed is well- supported and fully operational, thus providing an ideal 
environment for conducting research. In this section, a brief description of the CSI Program is 
first given, and then a description of the CEM test article and development of the analysis model 

follows. 

The NASA CSI Program has been formed with the overall objective of developing and 
validating the technology needed to design, verify and operate spacecraft in which the structure and 
control system interact beneficially to meet the requirements of 21st-century spacecraft [4]. 
Traditional spacecraft control for maneuvering is designed for a considerably stiff structure. 
Problems may arise if this same type of control is attempted to maneuver large flexible space 
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structures (LFSS). Interaction between the structure and control system can reduce performance, 
restrict operations, and cause dynamic instabilities. The interaction between controls and structure 
may also be detrimental when the flexible appendages such as solar arrays or antennas are 
articulated. Hence the structure and control interactions must be considered in order to achieve 
desired system performance. 

A CSI testbed has been developed in the Spacecraft Dynamics Research Lab (SDRL) at 
Langley Research Center. Some objectives of the testbed are to validate CSI design methodology, 
implement practical sensors and actuators for use in LFSS control, and to evaluate controller 
designs. This testbed, shown in Figure 2, consists of the CSI Evolutionary Model (CEM), 
suspension system, instrumentation for sensing, actuation, and data acquisition and processing 
capabilities. The CEM was designed to possess dynamic properties typical of space platforms 
proposed for remote sensing and communications. Line-of-sight (LOS) pointing is the primary 
objective for control system designs. 

The CEM consists of a 52.5 ft. long main truss, a 16 ft. diameter reflector, and a 9.2 ft. 
long laser tower as shown in Figure 2. To monitor the LOS pointing accuracy, a low powered 
laser has been mounted on the laser tower such that the beam reflects upon a mirrored surface 
mounted on the reflector. A photo-diode array above the reflector measures the beam reflection to 
give LOS accuracy of the CEM [5]. There are 62 bays in the central truss, and the tower truss has 
1 1 bays. Each bay is a 10 inch cube made of aluminum struts with single-laced diagonals. The 
reflector consists of 8 tapered ribs attached to a baseplate and bent to a preformed shape by a 
tension cable at their tip. The structure is suspended from the ceiling by two steel cables 64.5 ft. 
long. The total structure weighs 741 lbs. 

For control and excitation purposes, eight actuators and eight accelerometers are placed in 
pairs at four locations along the truss (Figure 3). The actuators are compressed air thrusters which 
deliver proportional bi-directional force. Each actuator is capable of producing 2.2 lbs. of force. 
The control sensors are servo accelerometers with 5 volts/g. sensitivity. 
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Figure 2: The CSI evolutionary model 
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The evaluation model used in this study is the most accurate mathematical model of the 
CEM currently available. For this reason it is also called truth model This model is used to predict 
the behavior of the CEM for various scenarios such as response to a disturbance or response to 
various closed-loop control algorithms. There are several steps involved in the development of the 
evolution model, and each step involves the development of a different type of model. First, a 
detailed finite element model is developed. This is also referred to as an analysis model. Next, 
data from system identification tests used to develop an experimental model. Finally, information 
from the experimental model is used to improve the finite element model. For large complicated 
structures such as the CEM, this cycle may be repeated a number a times. Each time the cycle is 
completed, the analysis model becomes a more accurate mathematical representation of the physical 
model. For this study, the analysis model is used in state-space form, which is further described 

later in this section. 

A. Finite element model 

A detailed finite element model (FEM) of the CEM (Figure 3) was developed using the 
NASTRAN finite element analysis program. The FEM is made up of beam, rod, and plate 
elements with 609 grid points and contains over 3000 degrees of freedom. The suspension cables 
are included in the FEM and are modeled by rod and spring elements. Lumped masses are also 
included to represent the joints and actuators. 

The analysis model predicts 8 1 modes of vibration below 50 Hz. Thirty of these modes are 
below 10 Hz [5]. Since the CEM is cable suspended, there are six rigid-body modes (also called 
pendulum modes) with frequencies between 0.1 Hz and 0.9 Hz. The first three flexible modes 
(modes 7, 8, and 9) involve bending and torsion of the CEM. These mode shapes are shown in 

Figure 4. 
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Figure 4: The first three flexible-body modes of vibration 


B. Experimental Model 

System identification tests are conducted to identify modal vibration parameters including 
frequencies, damping and mode shapes. These identified parameters define the experimental 
model. System identification for the CEM consists of multi-input, multi-output (MIMO) tests, and 
extraction of the frequency response functions (FRF) between the acceleration output to force 
input. Simply defined, an experimental FRF is found by dividing the Fast Fourier Transform 
(FFT) of the output by the FFT of the input. Typical FRFs from the center of the main truss, in 
the horizontal and vertical planes, are shown in Figure 5 [5]. Also shown on these plots is the 
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some of the dominant modes; however refinement of the FEM is needed The modal vibration 
parameters for nine modes have been identified at the present time. These parameters were used to 
refine the FEM and is included in the current evaluation model. Further system identification tests 
are underway to better identify the modal parameters for flexible modes below 10 Flz. 


Y - FREQUENCY RESPONSE FUNCTION 




C. State-space model 

The state-space form of the analysis model contains the first 40 modes from the refined 
FEM. The process of transforming the FEM to state-space form is now described. 

An FEM with n physical degrees-of-ffeedom (dof), m actuators (inputs) and m sensors 
(outputs) can be represented by the following second order dynamic system. 


Mz + Cz + Kz = Qu and y = Pz 


( 1 ) 
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where M is an n x n mass matrix, C is an n x n damping coefficient matrix, K is an n x n stiffness 
matrix, and z is an n x 1 vector containing the physical dof of the model. Further, u is an m x 1 
input vector, Q is the configuration matrix for the actuator locations, y is an m x 1 output vector (in 
terms of acceleration since the output sensors are accelerometers), and P is the configuration matrix 
for the sensor locations. 

Note n system eigenvalues A., and n corresponding eigenvectors d, are calculated using 
finite element analysis. The number of modes to be kept for the state-space model is / and is less 
than n. Thus the modal matrix V has as its columns / n x / eigenvectors, and is orthogonal with 
respect to M and K. Hence V can be normalized so as to satisfy 

V T MV = I , and V T KV = A ( 2 ) 

where A = diag[Xj] = diag [cof]. coj (i = 1,2,...,/) are natural frequencies of the system. If a vector 

is defined as z = Vq, where q is the l x / modal coordinate vector, then substitution by the modal 
vector into Eq. (1) and premultiplication by V T results in the modal system 

q + C'q + Aq = Q'u 

y = pvq = - PV(C'q + Aq) + PVQ'u ( 3b ) 


where C = V T CV is a real symmetric l x / matrix, and Q' = V T Q is a real / x m matrix. The system 
is now in reduced modal coordinates. 

If an / x / state vector x is defined as x = [q T q T ] T , Eq. (1) can be rewritten in the 
following state-space form: 


where 


x = Ax + Bu and y = Cx + Du 
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C = [- PVA -PVC'] and D = PVQ' 


(5b) 
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The state matrices A, B, C, and D are 21 x 21, 21 x m, m x 21 , and m x m, respectively. The state- 
space form of the analysis model is used extensively in this study. A further treatment of the 
transformation from a FEM in physical coordinates to a state-space modal model is found m [6]. 
The actual values of (A, B, C, D) for the CEM are given in [7]. 

IE. Selection of Filter Specifications 

Before a filter is designed, a set of critical frequencies and the attenuation desired at these 
frequencies must be specified. These filter specifications determine the filter order, and in some 
cases, the cutoff frequency (such as with Butterworth and Bessel filters). The selection of these 
specifications depends on the frequency characteristics of the signal to be filtered, and the 
application of the filtered signal. For the reference problem, low-pass filtering of accelerometer 

signals to be used in closed-loop control applications is desired. 

Frequency response functions (FRFs) from system identification tests are examined for 
frequency content of CEM accelerometer outputs. A large number of significant spikes appear in 
the 0 to 10 Hz range. Some of these represent the pendulum and structural modes to be controlled. 
There are a number of significant spikes, such as one at 34 Hz, as well as many smaller spikes 
above 10 Hz. Some of these are structural modes, and many are attributed to high frequency noise 

in the accelerometers. 

In preliminary closed-loop control tests, a 34 Hz mode became unstable in the attempt to 
control nine modes with frequencies below 2 Hz [10]. A possible solution for this problem is to 
use a low-pass filter to attenuate this 34 Hz frequency component This suggests a stopband 
frequency of 30 Hz. A significant portion of the frequency components above 30 Hz needs to be 
rejected so a 90% loss (20 dB) is specified as the minimum attenuation desired in the stopband. 

For this problem, the passband frequency and maximum attenuation are more difficult to 
specify. Ideally, all frequency components above the frequency of the highest controlled mode 
(1.9 Hz for the CEM) would be rejected. Thus, the passband frequency for an ideal filter is 2 Hz. 
However, since practical filters exhibit rolloff near the passband frequency, the passband 
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frequency needs to be significantly higher than 2 Hz. Also, the phase loss near the cutoff 
frequency could significantly distort the frequency components of the modes to be controlled. If 
the passband frequency is specified as 15 Hz, there is minimal phase loss in the region below 2 
Hz; also the phase response of the filter in this region is fairly Unear. Maximum attenuation in the 
passband is specified as 0.5 dB since accuracy of the frequency components of the controlled 
modes is desired. For the Butterworth and Bessel filters, the cutoff frequency is specified as 15 
Hz. Although there are 3 dB of attenuation at the cutoff frequency for these two filters, with the 
cutoff frequency set at 15 Hz, attenuation in the region below 2 Hz is less than the specified 0.5 
dB. Now the filter specifications can be designated in the foUowing manner: 

f p = 15 Hz , <x max = 0.5dB, f s = 30 Hz , ct mi „ = 20 dB 

and are iUustrated in Figure 6. The dashed curve represents the frequency response of a low-pass 
filter that meets the specifications. For the numerical simulations, Butterworth, Chebyshev, 



Figure 6: Filter specifications for CSI evolutionary model 
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Cauer, and Bessel filters are designed according to these specifications, and used to filter the 
simulated acceleration output from the CEM. 

IV. Open-Loop Simulation Results 

The n um erical open-loop tests (system plus filter) are conducted to extract the frequency 
response of the system when subject to an impulse force. The frequency response functions for 
the open-loop system are calculated for fifteen CEM and filter combinations. The first step is to 
compute the fiequency response of the system without any sensor filters. This response is 
compared to frequency responses of all the other system configurations. The other system 
configurations include the CEM with either a low-pass Butterworth, Chebyshev, Cauer, or Bessel 
filter. The frequency range examined is between 0 and 50 Hz. The resulting complex responses 
are plotted in phase and magnitude form. 

The first set of frequency response functions show the effect of using second, fourth and 
sixth-order Butterworth filters with the CEM. The magnitude and phase of these frequency 
responses, and that of the unfiltered system is plotted in Figure 7. The magnitude responses in the 
region below 10 Hz are indistinguishable from each other. The frequency requirements below 2 
Hz are satisfied by each filter. The magnitude response of the systems with fourth and sixth-order 
Butterworth filters fulfill the requirement that the magnitude of the frequency components above 30 
Hz be attenuated by 20 dB (90% loss). Thus a fourth-order Butterworth filter is needed to satisfy 
the filter specifications. The phase loss caused by the Butterworth filters increases significantly 
with filter order. In the region below 2 Hz the loss appears to be linear. 

With a system with many modes it is somewhat difficult to specify relative stability of a 
system using the phase and gain margins. However, trends in the frequency responses of the 
systems with various filters can be identified. There is good relative stability of the system without 
a filter, and even with the second-order Butterworth filter. This is illustrated by examination of 
Figure 7. The phase response in Figure 7b for the system without a filter, and with the second- 
order filter never exceeds 180 degrees in the frequency range shown. Thus there is some positive 


Phase (deg) Magnitude (dB) 


(a) Magnitude Response 



(b) Phase Response 



Figure 7: Frequency response 
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phase margin at any gain crossover frequency, and the gain margin is infinite. However, one can 
see that the addition of any filter, especially the higher-order filters, causes the phase margin to 
decrease significantly. The phase crossover frequency of the system with the fourth-order 
Butterworth filter is approximately 16 Hz, and that of the sixth-order filter decreases to 
approximately 12 Hz. The general trend is that the relative stability deteriorates with increased 

filter order. 

The second set of frequency response functions show the effect of using second, fourth 
and sixth-order Chebyshev filters with the CEM. The magnitude and phase of these frequency 
responses, and that of the unfiltered system is plotted in Figure 8. A third-order Chebyshev filter 
satisfies the specified requirements of at least 20 dB loss in the stopband, and a maximum 
attenuation of 3 dB in the passband. When the maximum attenuation is reduced to 0.5 dB, the 
slope of the magnitude response decreases in the transition band. Thus for a third-order 
Chebyshev filter with 0.5 dB maximum attenuation in the passband, the attenuation at the stopband 
frequency (30 Hz) is less than 20 dB (approximately 19 dB). To satisfy these specifications, a 
fourth-order Chebyshev filter is needed. The fourth-order Chebyshev filter meets the required 

specifications as seen in Fig. 8a. 

The trend of deteriorating relative stability of the responses of the systems with the fourth 
and sixth-order Chebyshev is similar to that of the Butterworth filters. The phase crossover 
frequency for the fourth-order Chebyshev system is approximately 14 Hz, and for the sixth-order 

system, 9 Hz. 

Figure 9 shows the frequency responses for the CEM with second, fourth and sixth-order 
Cauer filters. The second-order Cauer filter satisfies the specifications. As with the Chebyshev 
filter, the decrease in maximum attenuation in the passband from 3 dB to 0.5 dB causes the slope 
of the phase response to decrease. A third-order Cauer filter is required to satisfy the specification 
that a minimum of 20 dB attenuation be attained by the stopband frequency (30 Hz). Although the 
frequency response of the CEM with a third-order Cauer filter is not shown, it is obvious that at 30 
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Magnitude Response 



Frequency (Hz) 


(b) Phase Response 



Frequency (Hz) 

Figure 8: Frequency response of CSI evolutionary model with Chebyshev filters 
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(a) Magnitude Response 



(b) Phase Response 
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Hz, the response of the system with a second-order filter is not 20 dB below that of the system 
without a filter. 

The frequency response of the systems with the Cauer filters are very similar to that of the 
Chebyshev systems. However, the attenuation rate in the transition band of the Cauer systems are 
slightly higher than with the Chebyshev system. Also, due to the ’ripple’ in the stopband of the 
Cauer filter, the magnitude response in the stopband is slightly more distorted. The same 
observations given about relative distortion of the Chebyshev systems can be given for the Cauer 
systems. The phase crossover frequency for the system with the four and sixth-order Cauer filter 

is approximately 14 Hz and 9 Hz, respectively. 

The last set of frequency response functions are for open-loop systems consisting of the 
CEM and second, fourth and sixth-order Bessel filters. The magnitude and phase responses for 
these systems are shown in Figure 10. A tenth-order Bessel filter satisfies the specification of 3 
dB attenuation at the passband frequency but does not attain 20 dB loss at the stopband frequency. 
As with the Butterworth filters, in the range below 2 Hz, the Bessel filters cause less than 0.5 dB 
attenuation of the CEM output. The sixth-order filter provides approximately 17 dB of attenuation 
(86% loss) at the stopband frequency. Since higher order Bessel filters do not provide much 
additional attenuation at the stopband frequency, if 20 dB loss at 30 Hz is critical, then a Bessel 
filter is not appropriate. 

As with the other open-loop systems including filters, relative stability deteriorates with the 
addition of Bessel filters to the system. The phase crossover frequency for the fourth and sixth- 
order filters are approximately 22 Hz and 18 Hz, respectively. The Bessel filter causes the least 
amount of phase loss, and has the most linear phase response of all the filters considered. 

The difference in the frequency content of a filtered and unfiltered accelerometer signal 
from actual laboratory tests is shown in Figure 11. In the two tests, several dominant modes of the 
CEM were excited and then allowed to decay freely. In the first test, sensor filters were not used. 
In the second test, second-order Cauer filters were used to filter accelerometer output. FFTs were 
performed on the free decay portion of the output. The magnitude of the frequency response is 
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found by calculating the absolute value of each FFT function. There is a good amount of activity 
in the upper half of the magnitude response in Figure 1 la. This can be attributed to noise, 
disturbances and other unknown factors. In Figure 1 lb, the portion of the magnitude response 
above the cutoff frequency ( 15 Hz) is much smoother compared to that of Figure 1 1 a. The Cauer 
filter has successfully filtered out the unwanted frequency components of the signal. 

By examining the frequency responses of the filtered and unfiltered open-loop systems, the 
effect of a sensor filter on the phase and magnitude of the CEM output can be seen. As higher- 
order filters are added to the system, the relative stability of the system decreases. This will greatly 
affect the ability to use the system in closed-loop control applications. 

V. Control Strategy 

Two objectives for control of large flexible space structures (LFSS) are accurate line-of- 
sight pointing and vibration suppression. Classical control designs depend on the accurate 
prediction of the dynamics of the structure. However, with large multi-bodied, complex 
structures, accurate models are difficult to develop. The models of LFSS available for control 
design usually contain inaccurate model parameters (such as modal damping coefficients). The 
problem of unmodeled high-frequency flexible modes is another challenge of control design. 
Although a detailed model for the CEM has been developed and preliminary identification has 
increased the accuracy of the first nine modes, there is still enough model error so that classical 
control methods are ineffective, or even worse, cause system instabilities. For these reasons, 
robust control designs are used. 

The two control designs used in this study utilize the knowledge of the design model to 
obtain maximum performance, but do not depend on model accuracy to maintain stability. The 
design model is the model which is used to design the controller. The first method, static 
dissipative control [8], is a local velocity feedback method, and has been implemented by the CSI 
Analysis Methods Group. The second method, virtual passive control [9], is a second-order 
controller using direct acceleration feedback and has been implemented by the CSI Ground Test 


Group [10]. 
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VI. Closed-Loop Simulation 

Sensor filters have been successfully used in the open-loop experiments discussed in 
Section IV. By examining the FFT curves in Figure 1 1 one can see how effectively a second-order 
Butterworth filter can condition the accelerometer signal. However, when closed-loop control of a 
LFSS is desired, filtering a signal can have a significant effect on the stability of the system. Due 
to this reason, experiments involving control of LFSS are usually conducted without extensive 
prefiltering of the sensor signal. Some systems may include a first-order filter, however, a first- 
oider filter is not very effective. With the larger, more complex, multibody space structures 
envisioned for the future, it may be necessary to include sensor filters so that closed-loop control 

of the structure can be accomplished more effectively. 

Numerical and laboratory experiments were conducted to investigate the effect of sensor 
filtering on performance of feedback controllers. All four low-pass filters, Butterworth, 
Chebyshev, Cauer, and Bessel have been initially investigated while the Butterworth and Cauer 
filters were used in further investigations. The Butterworth filter is used because it is a popular, 
easily implemented filter, incurring relatively conservative phase distortion. The Cauer filter is the 
other low-pass filter selected for study because it provides the best magnitude response 
characteristics; therefore it would be ideal to use for our purposes if the resulting phase distortion 

could be managed. 

Two types of numerical analysis are conducted on the closed-loop systems. The first type 
of analysis consists of calculating the eigenvalues of the system. The second type of analysis 
consists of calculating the response of the system to a given excitation and examining the time 
histories of these responses. Presented in this section are the results of numerical analysis of the 
closed-loop systems. The results of laboratory experiments are discussed in Section VII. 

The experiment to be simulated, consists of excitation of the structure for seven seconds, 
free decay for three seconds, and then closed-loop control is initiated to damp out vibrations of the 
system. Discretized simulations are conducted in an attempt to produce results that will match that 
of laboratory experiments. A sampling rate of 250 Hz is used because this provides a folding 
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frequency of 125 Hz which is well above the frequency of the highest predicted structural mode 
(approximately 46 Hz). 

The excitation consists of sinusoidal force input at four locations, one at 1.7 Hz (bending 
mode), one at 1.9 Hz (bending/torsion mode) and two at 0.145 Hz (first pendulum mode). 

Control force is applied at the eight thruster locations. Acceleration output is calculated for all eight 
accelerometer locations. Simulated acceleration at sensor location 1 (Figure 3) is shown in all the 
time history plots. 

The plots show acceleration output using the Static Dissipative Controller (SD) and Virtual 
Passive Controller (AVA) for various filters. The notation in the captions denote the filter type and 
filter order. 

In Figure 12b, controller SDB4 is used to control a system with fourth-order Butterworth 
filters. The response from each closed-loop system is stable, including that the response from 
systems with the sixth-order filters (Figures 12c and 13c). However, the performance of the 
controllers in the systems with the higher-order filters are obviously reduced. This is especially 
true for the static dissipative controllers developed using the sixth-order Cauer filter (SDC6), and 
sixth-order Butterworth filter (SDB6). 

In Figures 15 and 16, the output of the system using controller AVA, and second and third- 
order Butterworth and Cauer filters are plotted. The unstable modes calculated with eigenvalue 
analysis (12.4 Hz, and 16.6 Hz) are only slightly noticeable in the systems with third-order filters. 
The eigenvalue representing an unstable mode is very close to the imaginary axis, so this mode 

could be considered only marginally unstable. 

In summary, all closed-loop cases illustrate that the filters perform this appropriate function 
and provide to the controller a clear signal to be controlled. The closed-loop program is stable and 
whatever phase distortion that is introduced by the filters have minimal effect on the controller s 
performance. This may not be true, in general, however, and an all-pass filter may be required to 
remove the nonlinear phase characteristics introduced by the filters. 
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VII. Conclusions and Future Activities 
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Numerical analysis shows that it is possible to use sensor filters in closed-loop control of 
large flexible space structures (LFSS). The dynamics of higher-order filters must be considered 
however because the filter phase rolloff causes phase distortion to the output signal. The stability 
margins of an open-loop system including filters decreases as the order of the filter increases. 

Active vibration suppression of a system including the CEM and sensor filters was 
demonstrated using numerical analyses. Two robust controllers were designed using the static 
dissipative method and the virtual passive method. Both were developed and evaluated for a 
system which included filter dynamics. The systems with the static dissipative controllers and 
second and third-order filters, and the systems with the virtual passive controller and second-order 
filters remained stable. The systems with higher-order filters became unstable when feedback 
control was attempted. When the filter dynamics were included in the design model for the static 
dissipative controller, the corresponding closed-loop systems remained stable. Numerical analysis 
of closed-loop systems with even sixth-order filters remained stable. However, performance of 
the controller was greatly decreased. 

Laboratory experiments were conducted to validate the numerical analysis. The controllers 
successfully decreased vibration of the CEM without sensor filters after excitation of several 
bending and pendulum modes. However, when second-order low-pass filters were implemented 
digitally to condition sensor output, a 7 Hz mode was excited and became unstable in the closed- 
loop experiment The static dissipative controllers which were developed with filter dynamics in 
the design model were implemented next and the same 7 Hz mode became unstable. Only the 
controllers developed with sixth-order filters were implemented successfully. These controllers 
have very small gains and therefore exerts only small amounts of force on the CEM. Although the 
higher-order mode was not excited, the controller did not provide much additional damping to the 

system. 

There are several reasons for the difference in results between the laboratory experiments 
and the numerical simulations, including model uncertainties and discretization errors. Further, 
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when the low-pass filters are combined with a physical system, the combined system is no longer a 
pure continuous-rime system, thus destroying some of the properties of the original system which 
were sufficient to guarantee stability of the closed-loop system. A sufficient condition for stability 
for both the static dissipative method and the virtual passive method is that the open-loop system's 
mass (M) and stiffness (K) matrices must be positive definite and the damping matrix (D) must be 
at least positive semi-definite for stability. When the filters are combined with the structure, it may 
not be possible to express the resulting system as a second-order system with a positive definite M 
and K, and a D that is at least positive semi-definite. This problem suggests one direction for 
further research: design a filter which can be combined with a physical structure such that the 
resulting system can be expressed as a stable second-order system. 

The problems with successful implementation of sensor filters for closed-loop experiments 
in the laboratory suggest several future activities. First, an investigation of how the digital filters 
are implemented may suggest changes in the real-time computer software. Also, further study of 
discretization methods for the filter state-space model needs to be made. The results of this study 
also suggest different approaches for designing a low-pass filter for closed-loop applications. 

One approach is to modify the transfer function of a conventional low-pass filter by adding 
zeros to the numerator in order to compensate for the phase lag caused by the filter. Another 
approach for low-pass filter design is to formulate the filter as a second-order system. This 
corresponds to the virtual passive approach for controller design. The filter would be designed as 
a virtual mass, spring and dashpot system. The mass, stiffness, and damping coefficients would 
be chosen so that when it is ’attached’ to the physical structure, the mass and stiffness matrices of 
the combined system would be positive definite, and the damping matrices would be at least 
positive semi-definite. 

Control law design for large flexible spacecraft is becoming more challenging as these 
structures increase in complexity. Typically such systems contain many flexible modes densely 
packed within a small frequency range. Also there are many unmodeled modes which may affect 
stability of the closed-loop system. Closed-loop control requires output from sensors for 



32 


information about the modes to be controlled Sensor filters which attenuate unwanted frequency 
components would greatly facilitate this process. The successful design of sensor filters which do 
not destabilize closed-loop systems may greatly simplify the control design process. 
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Appendix C 

Start-Up Procedures and Hardware/Software Notes 


I: Operational (Start-Up) Procedures. 


This system assumes the user has an MS-DOS environment using an IBM-compatible AT 
286-based computer (upgrades of 386 or 486-based machines would require minor revision of the 
software). Further the software requires the Franklin C-51 compiler (see software license 
agreement section below) and the Dallas software support for the DS2250 (see license agreement 
section below). Both of these packages are on the floppy disc provided to NASA Langley 
Spacecraft Dynamics Branch. 

Upon entering the DOS environment, simply load file ”tf.Bat" and execute. This file 
prompts the user for transfer function coefficients (stored in an include file). Next the control 
program is compiled and downloaded to the microcontroller. The load controller (DS2250) takes 
over and imm ediately execute the control program. 

II. Software Limitations 

Continuous-time transfer functions are limited to fifth order (due to the direct look-up table 
for the bilinear-z transformation) and discrete-time transfer functions are limited to 15th-order. To 
change the discrete-time transfer function order, simply modify the DEFINE statements located at 
the end of the tf.BAT file to include more coefficient terms (see Appendix A). To change the order 
of the transfer function in the continuous-time case is more involved. It may be easier to perform 
the discretization operation off-line (using MATLAB, for example), and enter the controller/filter in 

discrete-time form. 

TTT. Hardware Issues 

The bandwidth of the DPU is subject to the data acquisition system and control algorithm 
being tested. The data acquisition system has a limiting constraint at the A/D end of the system of 
lOKHz (the D/A is much faster). The DS2250 has a 12MHz clock. A +5v supply and a ±12v 
voltage reference power supply are required to power the DPU. 


A switch is located on the DPU to transfer control from the supervisor system to the load 
controller. Be sure that this switch is turned OFF when downloading and/or when SW1 is in the 
stop position (see Figure 11). 
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IV. Site License Agreements 


In order to fully employ the DPU software, site license agreements must be purchased. 
Please contact: 


(a) Franklin Software Customer Service 
888 Saratoga Avenue, Suite #2 
San Jose, CA 95129 

Phone (408) 296-8051 

fortheC-51 C-compiler 


(b) Dallas Semi Conductor 

4401 South Beltwood Parkway 
Dallas, Texas 75244-3292 

Phone (214) 450-0448 

for the DS2250 software or upgrades if desired. 


